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THEORY OF THOMSON SCATTERING FROM A WEAKLY IONIZED PLASMA* 

by Richard G. Seasholtz 
Lewis Research Center 

SUMMARY 

The theory of Thomson scattering from a weakly ionized plasma is extended to in- 
clude the effects of unequal electron and ion temperatures, a constant magnetic field, and 
electron drift. The Born scattering formula and the fluctuation-dissipation theorem are 
used to calculate the differential backscattering cross section. The plasma is described 
by the hydrodynamic (continuum) equations for the conservation of mass, momentum, and 
energy of each species. Included are electron-neutral and ion-neutral collisions but not 
collisions between charged particles. The results given are valid if the wavelength of the 
probing wave is greater than the ion-neutral mean free path. With no magnetic field the 
backscattered power is shown to be generally less than that found from a collisionless 
calculation and to be dependent on the ion-neutral and electron -neutral collision frequen- 
cies when the electron to ion temperature ratio T e /T^ does not equal 1. A magnetic 
field normal to the incident wave causes the backscattered power to be increased for 
T /T. in the common experimental range of 1 to 10; it is shown that the scattered power 
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may increase with T g /Tj rather than falling off in accord with the usual (1 + Tg/Tj) - 
relation. The effect of electron drift (less than the critical velocity corresponding to the 
onset of instability) on the Thomson scattering spectrum and on the total scattered power 
is calculated. As the electron drift velocity is increased, the spectrum develops a sharp 
peak at a frequency shift corresponding to the ion acoustic wave, and the total scattered 
power is enhanced. 


INTRODUCTION 

An important new technique of plasma diagnostics is based on Thomson, or incoher- 
ent, scattering of electromagnetic radiation. Although originally applied to the study of 
* 

Based, in part, on a dissertation submitted to Case Western Reserve University in 
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the ionosphere (refs. 1 to 6), Thomson scattering measurements have also been used in 
the diagnosis of laboratory plasmas using high power lasers as the radiation source 
(e.g. , ref. 7). The measurement of the Thomson scattered light of a pulsed ruby laser 
has been suggested (ref. 8) as a method of studying re-entry plasmas. Other possible 
applications of Thomson scattering of laser light include measurements in rocket exhaust 
and explosion generated plasmas. 

Thomson scattering measurements offer the important potential advantage of being 
able to measure many quantities of interest, such as electron and ion densities, temper- 
atures, drift velocities, and even velocity distributions at a specific interior point in the 
plasma. Techniques that involve the insertion of physical probes (e.g. , Langmuir 
probes) into the plasma offer much less spatial resolution and can disturb the plasma. 
Other radiation based diagnostic methods either measure what leaks out of the plasma, 
as with spectroscopic methods, or measure an average value along a path through the 
plasma as with microwave and laser interferometry techniques. 

Thomson scattering measurements were first used in studies of the ionosphere. In 
1958, Gordon (ref. 1) suggested that, if a beam of radio waves with a frequency well 
above the plasma cutoff frequency were sent into the ionosphere, the radiation scattered 
by the free electrons could be detected. Analysis of the spectrum and scattered power 
could provide information about the electron density and temperature. He noted that this 
measurement could be made above, as well as below, the altitude of maximum electron 
density. In his calculation he assumed that the electrons were in random thermal motion 
similar to neutral particles. He predicted that the scattered signal should have a spec- 
tral width corresponding to the thermal velocity of the electrons since their motion would 
cause the scattered signal to be Doppler shifted with respect to the incident wave. 

Shortly after Gordon's initial work, Bowles (ref. 2) successfully observed this phe- 
nomenon using high power radar. The scattered power agreed with that predicted by 
Gordon; however, Bowles was not able to observe the broadening of the spectrum that 
had been calculated. He explained this discrepancy by modifying the theory to include 
the electrostatic interaction between the electrons and ions and concluded that the spec- 
tral width should correspond to the ion, rather than to the electron thermal velocity. 

With the observations of Bowles as a stimulus many workers attacked the problem 
theoretically. The early work was done by Dougherty and Farley (ref. 3), Salpeter 
(refs. 9 and 10), and Fejer (ref. 11), who dealt primarily with plasmas in thermal equi- 
librium. More recent works by Salpeter (ref. 12), Rosenbluth and Rostoker (ref. 13), 
Buneman (ref. 14), and Farley (ref. 6) have extended the theory to include unequal elec- 
tron and ion temperatures. All this work has been based on kinetic theory calculations. 

Cohen (ref. 15) took a somewhat different approach and used very simple two-fluid 
collisionless hydrodynamic (continuum) equations to calculate the Thomson scattering 
spectrum. His continuum theory disagreed with the kinetic theory calculations in that 
the spectrum he obtained consisted of only sharp spikes since the broadening effect of 
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Landau damping was not included. A more complete continuum theory based on both two 
and three fluid models was given by Tanenbaum (ref. 16) for equal electron and ion tem- 
peratures. His calculations included collisions between the charged particles and the 
neutrals and gave results that agreed very well with the calculations of Dougherty and 
Farley (ref. 5) who had included collisions in their kinetic theory. 

The advantages of the continuum theories over the kinetic theories are that they are 
inherently simpler and collisional effects are easily included. Indeed, the continuum 
theories are only valid when the mean free path is less than the characteristic length of 
the probing wave and the collisional damping outweighs the Landau damping. 

In this report the continuum theory of Thomson scattering is extended to include 
scattering from plasmas with unequal electron and ion temperatures, a constant magnetic 
field, and electron drift. Electron-neutral and ion-neutral collisions are included in the 
theory but collisions between charged particles are neglected. Thus the theory is not ap- 
plicable to plasmas with very low electron temperatures where electron-ion collisions 
become dominant. The method used follows that of reference 3 (also see ref. 17), which 
used the Born scattering formula and the fluctuation -dissipation theorem which gives the 
spectrum of the fluctuations of a linear system in thermal equilibrium. The new results 
presented herein make it possible to analyze Thomson scattering experiments in the 
lower ionosphere and in other weakly ionized plasmas where collisional effects are im- 
portant and where kinetic theory calculations have not yet been made. Furthermore, in 
those cases where kinetic theory calculations have been made and where the continuum 
theory is also valid, the continuum theory has the advantage of being easier to apply. 


GENERAL THEORY 

The phenomenon of Thomson scattering in a plasma is dependent on fluctuations in 
the number density of the electrons. If the electrons were stationary and randomly dis- 
tributed, the total scattered power would be zero; however, because of the thermal 
motion of the electrons, fluctuations in the electron number density do exit. Further- 
more, when the Debye length is less than the wavelength of the probing wave, collective 
plasma effects come into play and the fluctuations of the ions influence the electron fluc- 
tuations. The scattering may be thought of as an interaction between the incident wave, 
with wave vector k , and a thermally excited wave in the plasma. The scattered wave, 
with wave vector k , is related to k by k = k + k (as shown in fig. 1) where k is 
the wave vector of the spatial Fourier component of the electron number density fluctua- 
tions that causes the scattering. 
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Figure 1. - Relation between wave vectors of incident, scattering, 
and scattered waves. 


Thomson Scattering Cross Section 

It can be shown (ref. 18) that the differential scattering cross section for a medium 
with a fluctuating dielectric constant is 

kj L 3 sin 2 p x 

°b( w o " ^ duJ = ^Ae(-k, co)j^dco (1) 

(4r r) 2 e 2 

where is the power scattered per unit solid angle per unit volume per unit frequency 
per unit incident power. (All symbols are defined in appendix A. ) In this formula co 
and k are the angular frequency and wave number of the incident wave, p is the angle 
between the electric vector of the incident wave and the direction of propagation of the 
scattered wave, co is the deviation of the frequency of the scattered wave from the inci- 
dent wave, and k is the wave vector of the scattering wave (1^= -2k o for backscattering). 
Also, ^|Ae(a)) is the power spectrum of the fluctuating part of the permittivity Ac and 
is defined as the Fourier transform of the auto correlation fluctuation of Ac; that is, 



(jAe(-k,o>)j 2 ) e +ic0T do) 


( A e*(-k,t)Ae + r)) 


(2) 


where the angular brackets on the right hand side denote a time average. Note that the 
k dependence of Ae arises because Ae has been expanded in a Fourier series in a box 
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of sides L, where L has been taken to be much larger than the scattering volume. 

To calculate Ae, recall that for transverse electromagnetic waves in a plasma with 
plasma frequency c o Q = (N Q e 2 /m e e 0 )^ 2 the dispersion relation may be written as 


,2 2 2 2 

k c = - u> 

o o e 


( 3 ) 


provided o> is much greater than the electron collision and cyclotron frequencies 

u 2 2 
(ref. 19). Using this dispersion relation for the plasma and noting that k^ = p Q e(r,t) gj^, 

the scattering cross section becomes 

a b (o> 0 - cd) dw = L 3 r 2 sin 2 p (| AN(-k, co)| 2 ) dw (4) 

2 2 

where AN is the fluctuating part of the electron number density and r 0 = e /4rr e Q m e c 
is the classical electron radius. Thus only the electron density fluctuations cause 
scattering. 


Fluctuation-Dissipation Theorem 

To use equation (4) the power spectrum of the fluctuations of the electron number 
density must be found. A form of the fluctuation-dissipation theorem due to Callen, 
Barasch, and Jackson (ref. 20), who generalized the original form of the theorem of 
Nyquist (ref. 21) to any linear system in thermal equilibrium is used. This theorem has 
previously been used by Farley, Dougherty, and Barron (ref. 4) to calculate Thomson 
scattering in a collisionless plasma. For a system in thermal equilibrium, the theorem 
may be stated as follows. 

Suppose that a set of generalized forces V^(t) are applied to the system. The system 
responses are represented by I^(t) where the spectral components of V. and I. are re- 
lated by 

v«) =S Y ij( w ) Vj(w) (5) 

and where Y.. is a generalized admittance tensor. The rate of change of energy of the 
system is 
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For high temperatures (or low frequencies) where kT » ho), the fluctuation-dissipation 
theorem may be stated as 

(lljM Ij(«)Ddc o = g (Y* + Y-.) do, (7) 

where the spectral function on the left hand side is defined by 



IjM I-j(co)l) 


e iwT 


do) 


= (l*(t) Ij(t + t)^ 


(8) 


that is, (fljM Ij (cc>) |) is the Fourier transform of the cross correlation function of the 
fluctuating quantities Ij(t) and Ij (t) . (For high frequencies where Kco > K T, eq. (7) must 
be modified to include quantum effects by replacing /<T with 1/2 K co coth (Ka'/2/<T) ; e.g. , 
for T = 300 K this change must be made if w> 4x10 rad/sec. ) 

The fluctuation-dissipation theorem is applied to a specific system which is in equi- 
librium by choosing the Ij to correspond to the desired macroscopic variable. The rela- 
tion for the rate of change of energy of the system is then used to identify the quantities 
that correspond to the V^. Finally, the equations of motion of the system are examined 
to find the proper admittance tensor that relates the I. and the V^. The fluctuation- 
dissipation theorem may then be applied to determine the fluctuations of the quantity cor- 
responding to the Ij. The fluctuations of the 1^ may also be thought of as arising from 
fluctuations of stochastic forces V . . The fluctuation-dissipation theorem may then be 
alternatively written as 

(|y.(o>) Vj(co)|) do> = “I (Y- 1 * + Yj: 1 ) dw (9) 

A simple example is to find the voltage fluctuations across a resistor R at tempera- 
ture T. The system is one dimensional, so the admittance tensor reduces to one element 
of value 1/R, and the fluctuation-dissipation theorem gives the power spectrum of the 
voltage fluctuations 


V(u>)| 2 ^ du> = — R dco (10) 

Note that since we allow co to take on both positive and negative values, this expression 
is 1/2 times the standard form of Nyquist’s theorem. 
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EFFECT OF COLLISIONS ON THOMSON SCATTERING WITH UNEQUAL 
ELECTRON AND ION TEMPERATURES (B = 0) 


The first case that is examined is the scattering of electromagnetic waves from elec- 
tron density fluctuations in a weakly ionized plasma with no magnetic field. Dougherty 
and Farley (ref. 5) have shown that for equal electron and ion temperatures, collisions do 
not affect the total backscattered power. Here it is shown that when the electrons and 
ions are not at the same temperature, the total backscattered power is dependent on the 
ion-neutral and electron-neutral collision frequencies. 

Farley (ref. 6) has calculated the scattering cross section for unequal electron and 
ion temperatures, but only for the collisionless case. Here it is found that for 
T e /Tj > 10, the total scattering cross section is substantially less than that of the colli- 
sionless calculation, but always greater than a cross section that varies with T /T. as 

_ 1 U 1 

[1 + (Tg/Tj)] . As the collision frequencies v en and zx n increase, the total scattered 
power approaches the (1 + T e /T^) - * curve, although it does not reach it even in the limit 
v Qn , i/j — Also, for T e /T^>5, the ratio of the collision frequencies v er J v \ n i 
i^ n fixed, has a large effect on the total scattered power, although the shape of the 
spectrum is essentially unchanged. 


Evaluation of Backscattering Cross Section With 
Unequal Electron and Ion Temperatures 

The backscattering cross section, from equation (4), is 

a b (u> 0 - w) da> = L 3 r 2 (|AN(2k Q , w) | 2 ) dco (11) 

so the scattering may be thought of as being an interaction between the incident wave, 
which is taken to be proportional to and the spatial Fourier component of the 

fluctuating part of the electron number density AN whose wave vector is 2k Q . (In the 
remainder of this report only backscattering is treated; however, the theory may easily 
be extended to other scattering angles by using equation (4) in place of equation (11) with 
k,= k g - k Q . The magnitude of k is then given by jkj = 2k Q sin e/2 where Q is the 
angle between k Q and k g . ) 

To calculate the thermal fluctuations of the electrons, a weakly ionized plasma is 
considered in which the electrons and ions each have Maxwellian distributions but not nec- 
essarily at the same temperature. It is assumed that the temperature of the neutrals and 
ions are equal. Because of the large difference between the electron and neutral masses, 
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the electrons tend to exchange energy slowly with the heavy particles so that each species 
may achieve a Maxwellian distribution, but at different temperatures. In order that this 
temperature difference may be sustained under steady-state conditions, it is hypothesized 
that some external energy source, such as a constant electric field, is acting on the 
plasma. 

The model used for the plasma is that of two fluids of charged particles, the electrons 
and ions, each of which is coupled to the neutral particles by collisions and which interact 
because of the self-consistent electric field. Short range electron-ion collisions are 
neglected. 

The fluctuations of the electrons and ions are calculated separately by applying the 
fluctuation-dissipation theorem to fictitious sets of particles that are identical to the real 
electrons and ions except that no forces due to the electric field are included. The actual 
fluctuations on the electron number density are then found by using Maxwell’s equations 
to include the interaction between the electrons and ions due to the self-consistent elec- 
tric field. 

In order to calculate the response function of each fictitious fluid, a longitudinal oscil- 
lating force F is applied to the s species (s takes on the value e for electrons and 
i for ions). The rate of change of energy of the s species in a volume L is 

— = L 3 F a • (12) 

dt ~s s v ' 

where = N u is the flux density of the s species. With no applied magnetic field, 
s s^s 

the response r g is in the same direction as the applied force F g , so the vector notation 
can be dropped. A scalar admittance function Y is introduced that relates the response 
to the applied force: 

r s = Y s F s < I3 > 

(Application of the fluctuation -dissipation theorem to a system with a tensor response 
function is treated in appendix B. ) These admittance functions will be obtained from the 
equations of motion of each set of fictitious particles. 

To apply the fluctuation-dissipation theorem, identify F with V and then, com- 
paring equations (6) and (12), note that the response I must be identified with L r and 
Y with Y L . The fluctuation-dissipation theorem (eqs. (7) and (9)) may thus be written 
in either of the forms 

(|r s | 2 }dco = — Re Y g dw (14) 

L 3 77 
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da) 


(15) 


The quantities F g and F^ may be considered to be stochastic forces that act on the 
electrons and ions, respectively. 

The self-consistent electric field, given by Maxwell's equations, will be an additional 
force that acts on the charged particles. The Maxwell equation 


VXH = !2 + J (16) 

at 

is used. Here V x H = 0, since only longitudinal oscillations cause fluctuations in the 
electron number density. 

Equation (16) may be expressed in terms of the flux densities as 


r e " r i = Yj(eE) (17) 

where 

Y t = iw — (18) 

a 2 

With the inclusion of the force due to the electric field, the equations for the electron and 
ions become 


r e = Y e< F e “ eE > 
r 4 = Y.(F. + eE) 


(19) 


These two equations together with equation (17) form a closed set of equations for the 
electron and ion flux densities that can be solved once we know F g , F., Y g , and Y^. 
Solving equations (17) and (19) for r g gives 


r e = 


F e Y e< Y i + Y I> + F i Y e Y i 
Y e +Y i + Y I 


( 20 ) 


Assuming that the stochastic forces on the electrons and ions are independent, that 
is, that ^F g (w) F^(cu) = 0, the power spectrum of r g is 
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(|r e (a>) \ 2 y dw = 


Y e< Y i + Y I> 
Y e +Y i + Y I 


(l F e(“)l 2 ) 


Y i 

Y i +Y I 


(|F i (w)| 2 )j dco (21) 


Using equation (15) for (|F e | 2 ^ and and the continuity equation for the electrons 

BAN 


Bt 


•+v.£ e = 0 


(22) 


yields the power spectrum of the electron number density fluctuations 


(]AN(k,co)| 2 )do> = — ^ 


2 T 3 

O) ttL 


Y e ( Y i + Y l) 


Y e + Y i + Y I| 



Y i 


i Y i +Y lI 2 


dco (23) 


To this point the analysis is essentially that of Farley (ref. 6) who then calculated the 
admittance functions for a collisionless plasma. Here Farley's work is extended by cal- 
culating the Y e and when collisions are important. For this purpose the approach 
of Tanenbaum (ref. 16) is followed and use is made of the transport equations for mass, 
momentum, and energy with the heat flow term and pressure tensor as derived by 
Goldman and Sirovich (ref. 22). Applying these equations to small amplitude disturbances 
in a weakly ionized gas with neutral velocity u n = 0 gives (to first order in the fluctuating 
quantities) 


9p s 

— - + P_V • U = 0 
Bt 30 ~ s 


(24a) 


3u 

~s + 

I^L- 

~s 

Bt 

Pso 

m 


s 

^s _ 

5 Pso 

0 p s 

Bt 

2 Pso 

Bt 


V 2 u +1 V (V • u_) 


1 os 


= - v u 
sn~ s 


J *aO 


2 X V 2 T = -N K V (T - T ) 

3 os s o K sn v s so ; 


(24b) 


(24c) 


The energy equations imply that there is some external force that causes the electron and 
ion temperatures to be different under steady -state conditions. The temperatures of the 
ions and neutrals are assumed to be equal. Note that s takes values e (for electrons) 
and i (for ions) and that n denotes neutrals. Also, p , q m p and T are the 
mass density, charge, mass, pressure, and temperature, respectively, for species s; 
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F is the force whose response we are calculating; p and p' = 2m p / (m + m ) 
are, respectively, the effective collision frequencies for momentum and energy transfer 
with the neutrals; and ij and \ , the viscosity and thermal conductivity for charged 

OS OS 

particles in a weakly ionized gas, are given by 


^os 


D p 
s sn 


^os 


15 *P S 

4m s C s^sn 


(25a) 


(25b) 


where D and C_ are constants (of order 1 or 2) that depend on the interparticle force 
s s 

law and the mass ratio between the electrons or ions and neutrals. 

It is assumed that p g = p go + p g , p g = p gQ + p g , and T g - T gQ + T g with the fluctua- 
tion quantities p', p', and T' proportional to eH^t-kz) Since longitudinal disturbances 
s s s 

have been specified, F g and u g are in the z direction. 

Using the mass and energy equations allows the momentum transport equations to be 
written in the following form: 


where 


iu,N o F s 
N u : 2. _! 

° s 2 z 
k\T g z s 


(26) 


Z s ~ A S + 2 l 0 S^S + 2 < 3 W 1 ] " 2e .q 


(27) 


1 + i (*a? 

3<t ) 


A s = 


1 + if"-’ 


(28a) 


5k 2 *T 


a s = v sn + 


2C s^sn m s 


(28b) 


and 0 and \p are the normalized Doppler shift frequencies and collision frequencies 
s s 
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(29a) 


*«,= 


(29b) 


Comparison of equations (26) and (13) shows that the admittance functions for the 
electrons and ions considered separately are 


and 


i " N o 1 
k 2 *T e z e 


icoN 


Y. = 

l 


1 


k 2 Kl\ z i 


(30) 


Thus the power spectrum of the electron density fluctuations (eq. (23)) may be expressed 
in terms of the macroscopic parameters of the plasma and the differential backs cattering 
cross section (eq. (11)) may be written as 


CT b (w o 


- o>) dco = N Q r 2 


2 

2 


[ia + z^ 

Im 

2 z. 

a x 

z + - 

« 2 (Zi + M z e ) + z e z i 

e 0 

^ l±a + z. 


dcr 

tto; 


(31) 


where p = T g /T^ and a = (kX^) - *, with the electron Debye length, X^, equal to 
( K ^e e o/^o e ^l/2' This is essentially the same expression, with slightly different 
notation as Farley's (ref. 6). However, the z g and z^ used here have been derived 
from continuum equations rather than from the kinetic equations Farley used. 

Examination of equations (27) and (29a) shows that z * (co) = z (-w). Any function 
F(z z.) with constant real coefficients also has this property, that is. F (cc) = F(-co)- 

6 7 1 / ey 

It follows that | F(o>) | z is an even function of w and that ImF(w) is an odd function of u\ 
Hence the differential scattering cross section (eq. (31)) is an even function of ce. 

For the case of equal electron and ion temperatures , the cross section simplifies to 
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(32) 


°bK 


gj) dco = -N r Im 
7 o e 


o ' 2 + z i 


a ( z i + z e> + z e z i 


dw 

7T0J 


Provided that the ion -neutral mean free path is less than the incident wavelength, equa- 
tion (32) gives results (Tanenbaum, ref. 16) that agree very well with those obtained by 
Dougherty and Farley (ref. 5) using the Boltzmann equation with the Bhatnagar -Gross - 
Krook (BGK) model for collisions with neutrals. 


Numerical Results 


The normalized differential cross section, as obtained from equation (31), is eval- 
uated for four values of (T e /T^) in figure 2 with a = 12. 7 and m. = 31 amu (this value of 
m. is chosen to represent a mixture of NO + and ions). Note that a typical Thomson 

scatter spectrum (ref. 3) has a peak at u> = kU (where U is the plasma sound speed) 

P P 



Normalized Doppler shift frequency, 6; = UWk)(mj/2/cTj) 1/2 

Figure 2. -Thomson scattering cross section as function of 
Doppler shift frequency, a = 12.7; normalized ion -neutral 
collision frequency, 1; normalized electron-neutral 
collision frequency, ip e , 0.1; ion mass, m } , 31 atomic mass 
units; Cj = Dj= 2; C e = D e = 1. 
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which is due to scattering from the so-called ion acoustic waves in the plasma. In addi- 
tion, the peak is sharp when there is little damping of the ion acoustic waves, but 
broadens and eventually disappears when the damping increases. Increasing the electron 
to ion temperature ratio causes the plasma sound speed to increase and the damping of 
ion acoustic waves due to collisions to decrease. This broadens the spectrum and in- 
creases the relative amplitude of the resonance, as can be seen in figure 2. The shape 
of the spectrum is essentially unaffected by a change in although the total power 

is changed. This is shown in figure 3 where the spectrum (normalized so that ct^(w 0 ) = 1) 
is plotted for two values of 


Electron to ion 
collision 



Figure 3. - Comparison of spectral shape for two values of normalized electron to ion 
collision frequency ratio. a= 12.7; normalized ion-neutral collision frequency, 
il/ j, 1; electron to ion temperature ratio, T e /Tj, 10; ion mass, m if 31 atomic mass 
units; Cj = Df = 2; C e = D e = L 


Discussion of Total Power 


The total backscattered power is given by 

°tot = f °bK “ w ) dco (34 

Although this integral may easily be calculated numerically (see fig. 5). an approximate 
analytic form can be derived that is valid for a » 1, ju < 10, and ^ > 1. Equation 
(34) may be written using equation (31) as 
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The first integral in equation (35) is easily evaluated using residue theory to find the 

Cauchy principal value, which is just given by -ni x (Residue at co = 0) since there are no 

poles in the lower half plane. Since z a ss 1, l a? + z l 2 may be approximated by 
o p e e 

(a + 1) , and the second integral can be evaluated in the same manner as the first. 

Thus equation (35a) becomes 
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If z = 1, the remaining integral gives zero contribution, and the total cross section 
is given by the first two terms in equation (36). This is equivalent to Farley's approxi- 
mate solution (ref. 6) and gives the identical result. This expression is also given by 
Salpeter (ref. 12) and Buneman (ref. 14). Note that there is no dependence on the colli- 
sion frequencies in this approximation. 

The case of the limit a — °° (Debye length —0) may be easily evaluated if it is as- 
sumed that A may be approximated by its low-frequency limit 

b 


A 


1 + 


2ico 



(37) 


Equations (28) show that the approximation is certainly valid for A g and should be fairly 
good for Aj when ^ > 1, provided p is not too large, since the main contribution to 
the total power then comes from the region < 1 (see fig. 2). Using equation (37) al- 
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lows the integral in equation (36) to be evaluated 


g tot _ l | U - i\ 1 

N j,2 1 + ii + 1/ I + Q 

o e 

where 


(38) 


with 


and 


Q 





(39) 



3I W 



K, 


1 + 


3D e^e 




-1 


(40a) 


(40b) 



(40c) 


For small Debye length ( a —°°), equation (38) shows that the total backscattered power is 
generally less than that calculated from the collisionless theory, but it is always greater 
than (1 + p)“* and is dependent on both the electron and ion collision frequencies. This 
approximate solution for x/s. = 1 is shown in figure 4 along with (1 + ju) - *, the numerical 
integration of equation (31) for a » 1, and the large T /T. approximate solution of 
Farley (ref. 6) for the collisionless case. For xf/. > 1 the approximate solution should 
be even closer to the exact solution. 

Figure 5 shows the total backscattering cross section (obtained by numerical integra- 
tion) for three values of i//. and two values of For large values of fj. the total 

backscattered power is less than that found from the collisionless theory, and it more 
closely follows the (1 + curve which is commonly used by experimentalists. The 

total power approaches the (1 + ju)“^ curve closely as becomes larger, although it 
does not equal it in the limit as and i//. — °° . 

It should be pointed out that \p e ayd !/a are independent of temperature only when 
the interparticle force is that of hard spheres (see ref. 23, p. 251). In general, this is 
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Figure 4. - Comparison of approximate solution (eq. (38}) and 
the true value obtained from the numerical integration of 
equation (31} for the total back scattered power. (Also shown 
are the approximate solution of Farley (ref. 6) for the collision- 
less case and the (1 + T e /T|) -i curve.) a = normalized 
ion-neutral collision frequency, ijj- v 1; normalized electron- 
neutral collision frequency, (jJ e , 0.1; ion mass, m,-, 31 atomic 
mass units; Cj = Dj = 2; C e = D e = 1. 


= (v sn /kKm s f 2 kT s ) 1/2 



Figure 5. - Total backscattered power as function of electron to 
ion temperature ratio from numerical integration of equa- 
tion (31). a= 12.7; ion mass, m } , 31 atomic mass units; 

C| ■ t>i * 2; c e = D e ■ 1 . 


not the case, and an additional temperature dependence would be introduced if the proper 
interaction were known and taken into account. However, this additional temperature 
dependence is not important to the calculations presented in this report. 


EFFECT OF COLLISIONS ON THOMSON SCATTERING IN A MAGNETIC FIELD 

In this section the theory of Thomson scattering is extended to include the effects of 
electron-neutral and ion-neutral collisions on the signal scattered from a plasma in a 
uniform constant magnetic field. 

It has been shown (refs. 4 and 10) that a magnetic field appreciably influences the 
Thomson backscattered signal only when the field is normal or nearly normal to the di- 
rection of propagation of the probing wave. Farley (ref. 6) has shown that for a collision- 
less plasma the total backscattered power is markedly altered by a magnetic field when 
the electron temperature T g is not equal to the ion temperature Tj. Since, for the case 
of no magnetic field, the scattered power has been shown to be dependent on the electron- 
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neutral and ion-neutral collision frequencies for T g 4 - T\ . it is not surprising that in the 
calculations presented here the total backscattered power depends on the electron and ion 
temperatures, the frequency of collisions, and also the magnitude and relative angle of 
the magnetic field. 

In the experimentally important range of T g /T^ between 1 and 10, the total back- 
scattered power is increased by the magnetic field - an effect also shown in Farley’s col- 
lisionless calculations. However, increasing the collision frequencies is shown to reduce 
the influence of the magnetic field and thus causes the scattered power to be closer to the 
field free case. It should be noted that when magnetic -field effects are important, the 
total scattering cross section is not proportional to (1 + Tg/T^ -1 as it is (to a good ap- 
proximation for small T e /T\) in the field-free case. 

Evaluation of Backscattering Cross Section With Magnetic Field 

Providing that the incident frequency is much larger than the electron plasma, colli- 
sion, and cyclotron frequencies, the differential backscattering cross section can then 
be written as (derived in appendix B) 


(o> - «) dw = — - Tz -1 • (B - B f ) • dw 

0 2 ria» L* * “ ~ J zz 


(B14) 


The coordinate system is chosen so that the incident wave vector k Q is alined with the 
z axis and the magnetic field B Q lies in the x-z plane and forms an angle (3 with k ( 
The symbol t denotes the Hermetian conjugate, and [ ] means that the expression is 
the zz component of the tensor . The tensors Z and B are defined by 


g = ge + ("If 1 + If 1 )" 1 

B = Z -iL^L/qZ: 1 + Z" i y i 

Po \ x x / 

11 V 1 


(B15) 

(B16) 


where a = T /T., Z and Z. are normalized impedance tensors that characterize the 
response of the electrons and ions, respectively, to an applied force proportional to 
e io>t-ikz^ and 2 ' 1 i s a tensor (closely related to the usual conductivity tensor) given by 
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where n = kc/cc. Note that, for equal electron and ion temperatures, B = Z and the 

rvi fvJ 

cross section (eq. (B14) is then just proportional to the imaginary part of the zz com- 
ponent of Z~*. 

It is only necessary to calculate Z„ and Z. (evaluated at k = 2kl to find the scat- 
he o 

tering cross section given by equation (B14). Other authors (refs. 4, 10, and 13) have 
used kinetic theory to calculate Thomson scattering when magnetic field effects are im- 
portant. In appendix C the continuum equations for a weakly ionized plasma in a constant 
magnetic field are used to calculate Z„ and Z. for the case where electron-neutral and 

ft; G £5^1 

ion-neutral collisions are important. 

An approximation used by other authors (refs. 4 and 10) may be used to simplify the 
amount of matrix algebra needed to evaluate the cross section as given by equation (B14). 
It is to assume that the phase velocity of any wave which causes scattering is much less 
than the velocity of light c. To use the approximation, the limit is taken as c ap- 
proaches infinity in the evaluation of equation (B14). The result of this is to decouple the 
transverse and longitudinal waves. (The longitudinal oscillations are the only ones that 
are important here since the scattering is caused by the fluctuations in the electron num- 
ber density. ) Applying this approximation leads to the following form of the backscatter- 
ing cross section 


ct. (w - to) do = N r'i |(Z _1 i z Im J 
b v o ' o e ' v a- 'zz 1 
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where 




-1 

zz 


a 


1 + 


p a 


Det(Z ) 

X 


& 


XX yy 


- Z 1 Z J 


xy yx' 


(43) 



z e z e - z e z e 

xx yy xy yx 
D et(Z) 


(44) 


and the elements of Z are 
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Z.. = zf. + faZT 1 + z; 1 ) 1 6- 6- 

i] i] i 'zz °1Z°]Z 


( 45 ) 


This is equivalent to the cross section used by Farley (ref. 6) for his kinetic theory 
calculations . 


Cross Section for p = 90° 


As discussed previously, the magnetic field affects the backscattering cross section 
only when the incident wave is normal or nearly normal to the magnetic field. In this 
section the special case of p - 90° is examined because at this angle the effect of the 
magnetic field is maximized. For other angles the full equations presented in the previ- 
ous section can be used if an evaluation of the cross section is desired. 

For p = 90°, the tensors Z defined by equation (C15) simplify to 





(46) 


where the upper sign is to be taken for the electrons and the lower sign for the ions . The 
quantities x . y z and A are now defined as 
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The thermal conductivity A g , given by equation (C14) becomes A og /[l + (0 S /^ S )^] 
for this special case of (3 = 90°. Note that the effect of a large magnetic field is to re- 
duce the thermal conductivity in the direction normal to the field lines. This change in 
the definition of a g (and hence A g ) may be made by replacing the C , which appears in 
equation (28b) with C [1 + (<p /rfs ) 2 ]. 


s L ' ^s' T s' 

Equation (46) may now be used to find Z and 


Z-, and these in turn may be used to 


evaluate the exact expression for the backscattering cross section by using equations 
(B14) to (B16). As in the general case, a simpler expression for the cross section may 
be found by using the approximation that the velocity of light is large compared with the 
phase velocity of the plasma waves that cause the scattering. With this approximation, 
the backscattering cross section for /3 = 90° becomes 
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where 
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nud y z , and A g are given by equations (48), (49), and (50). 

When the magnetic field approaches 0, A . = 0 and P = 1, thus equation (51) reduces 
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to equation (31) developed for the cross section with no magnetic field. Furthermore, for 
/3 = 0° (k | ) B q ), P g = 1 and hence the magnetic field does not change the scattering cross 
section when the magnetic -field lines are parallel to the incident wave vector. 


Spectral Magnitude at Zero Doppler Shift for p= 90° 

In order to better understand the effect of the magnetic field on the Thomson scatter- 
ing spectrum, the magnitude of the spectrum can be evaluated at zero Doppler shift for 
/3 = 90°. The case considered is that where the cyclotron radius of the electrons is much 
smaller than both the electron-neutral mean free path and the wavelength of the incident 
wave X Q but where the ion cyclotron radius is much larger than both the ion -neutral mean 
free path and X . In terms of the normalized variables defined by equations (C6), (C7), 
and (C8) these conditions may be expressed by 

4 » 1 

and 

4 « 1 


Furthermore, it is assumed that the electron and ion temperatures are equal and also that 

9 

a » 1. Then equation (51) is evaluated in the limit as co — 0 to find 
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Several observations can be made about the spectrum based on this expression. First, 
with ~ 0. 1 and ip- ~ 1, 0y,(cO is primarily a function of \p and <f> (or <f>. 

since <±> = (m./m ) L ' <b-) if <b d>, > 1. Increasing the magnetic field strength causes 
a b (cu o ) to become larger, which implies a narrowing of the spectrum since the total cross 
section is independent of both collisions and the magnetic field for T g = Tj. With zero 
magnetic field, increasing the collision frequencies causes the spectrum to be narrowed; 
however, here, where magnetic field effects are important, increasing the collision fre- 
quencies (up to the point where the ion collisions begin to dominate) causes the spectrum 
to be broadened. The physical explanation is that, as the electron-neutral collision fre- 
quency is increased, the electrons are somewhat freer to move across the field lines 
since their spiraling motion is broken up by the collisions; however, as the collision fre- 
quency is further increased, the collisions eventually dominate the effect of the field, and 
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the spectrum again becomes sharply peaked (as when B q = 0). Also the viscosity effects 
on the electron motion disappear for this case. The term involving the thermal conduc- 
tivity remains, although the role of xj/ is replaced by <p . Thus, the continuum equa- 

“ “ 2 2 
tions for the electrons retain their validity even for x^ e < 1, provided that <fi e + > 1. 

Finally, note that o^(o; o ) is a sensitive function of ^ e A^- This differs from the field- 

free case (for oft » 1) where has only a slight effect on the spectrum (refs. 5 

and 16). 


Numerical Evaluation of Cross Section 


In figure 6 equation (51) has been numerically evaluated to show the influence of col- 
lisions on the spectrum for (3 = 90°. The values chosen for the other parameters are 
typical of conditions found in the ionosphere at altitudes of about 100 kilometers if the in- 



figure 6. - Normalized Thomson scattering cross section 
as function of Doppler shift frequency for various values 
of normalized ion-neutral collision frequencies 
{$■ =(Vj /k)(m;/2/cTj) 1/2 ). Normalized ion cyclotron fre- 
quency ipj = ieB 0 /kmjHmj/2*cTj) 1/2 , 0.1; angle (£) between 
the incident wave vector and the magnetic field, 90°; a * 
12.7; electron to ion temperature ratio, T e /Tj, 1; normal- 
ized electron -neutral collision frequency, ip e , 0. 1 ipy, 
ion mass, m:, 31 atomic mass units; Cj * Dj = 2; 



cident wavelength A Q equals 6 meters. As discussed, we see that as i/a increases, the 
spectrum is broadened up to about )/a = 5. Then as i/a is increased further, the collisions 
begin to dominate the magnetic field, and the spectrum begins to peak about the origin. 

The effect of the angle /3 on the spectrum is shown in figure 7 for two different col- 
lision frequencies. This shows that the magnetic field affects the spectrum only for /3 
within a few degrees of 90° and that at the higher collision frequencies the relative change 
due to the field is less. 



Figure 7. - Normalized Thomson scattering cross section as function of Doppler shift frequency for various values of (3. Other 
parameters are the same as those in figure 6. 


For equal electron and ion temperatures , neither collisions nor magnetic field affects 
the total cross section, which for this case is given by reference 4 as 
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= Nr 
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o e 
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However, when the electron and ion temperatures are not equal, both collisions and 
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magnetic field strength influence the total cross section. In figure 8 the effect of a mag- 
netic field normal to the incident wave on the total cross section calculated by numerically 
integrating equation (51) is shown. This figure points out that in the relatively common 
experimental range when T /Tj is between 1 and 10, the total cross section increases 
with B q . Also, for large B field (electron cyclotron radius on the order of, or smaller 
than, the wavelength of the incident wave), the total cross section may increase with 
T /T. rather than falling off in accord with the usual (1 + T _/T-)~* relation. Farley 
(ref. 6) has published curves for the case of no collisions which show a similar effect of 
the magnetic field. 

In figure 9 the effect of collisions on the total cross section in the presence of a mag- 
netic field (0. = 0. 1) is shown. For Tg/T^ between 1 and 10, the effect of collisions is 
to drive the curve of total scattered power toward the field-free case. 

The effect of the angle fi on the total cross section is shown in figure 10. This 
again points out the fact that the magnetic field has an appreciable effect on the cross sec- 
tion only when the incident wave is nearly normal (/ 3 > 80°) to the magnetic field for 
T e /T t < 100. 


Normalized ion 
cyclotron 
frequency, 



Figure 8. -Total back scattered power as function of elec- 
tron to ion temperature ratio for eight values of nor- 
malized ion cyclotron frequency, normalized ion- 
neutral collision frequency, tjj j, 1 ; other parameters 
are the same as in figure 6. 


Ion cyclotron ..ngle 
frequency, 3 


0. 1 90 ° 



Figure 9 . - Total bacicscattered power as function of elec- 
tron to ion temperature ratio for several values of 
normalized ion-neutral collision frequency. Other 
parameters are the same as those in figure 6. 
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Figure 10. - Total backscatterecS power as function of 
electron to ion temperature ratio for various values of 
the angle between the incident wave and the magnetic 
field. Normalized ion cyclotron frequency, <£j, 0. 1; 
normalized ion collision frequency, 0j, 1. Other 
parameters are the same as those in figure 6. 

SCATTERING FROM A PLASMA WITH ELECTRON DRIFT 

Thomson scattering from a plasma in which the electrons have a net drift velocity 
relative to the ions has been discussed by Rosenbluth and Rostoker (ref. 13) for a colli- 
sionless plasma. They calculated the scattering for relative drift velocities below the 
critical value at which the instability of the ion acoustic wave occurs. They found that, 
as the drift velocity approaches the critical velocity at which the instability occurs, the 
spectrum develops a peak that becomes infinite as the critical velocity is reached. (The 
appearance of this infinity in the spectrum results from a breakdown of the linear theory. ) 

Here, the work of Rosenbluth and Rostoker is extended to include the effect of colli- 
sions, and backs cattering is calculated (in the stable regime) using the linearized contin- 
uum equations and the fluctuation -dissipation theorem. It is assumed that the electrons 
and ions have Maxwellian velocity distributions and, for the sake of simplicity, that the 
ions have zero net velocity. Our model is thus a plasma in which the ions and neutral par- 
ticles have no net fluid velocity and the electrons, which are assumed to have a Maxwellian 
velocity distribution, are drifting in the direction of the incident wave with a velocity be- 
low the critical velocity at which the two-stream instability occurs. Also, it is assumed 
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that the thermal velocity of the electrons is much greater than the drift velocity. In gen- 
eral this means that, for drift velocities close to the critical velocity, the electron tem- 
perature must be greater than the ion temperature. This seems to be reasonable, since 
for many physical situations in which the electrons are drifting with respect to the ions 
(e. g. , in a gas discharge device) the electrons will be hotter than the ions. 

The critical drift velocity is dependent on the electron to ion temperature ratio, the 
collision frequencies, and the magnitude and relative angle of the magnetic field. In gen- 
eral, the critical velocity is found to be relatively independent of the electron temperature. 
With higher collision frequencies (which increase the net damping) larger drift velocities 
are required to produce instability. If a magnetic field is present, the critical drift veloc- 
ity is smaller if the electrons are drifting perpendicular to the field. 

When the wavelength is much smaller than the Debye length, (or « 1), the only 
change in the spectrum of the scattered signal is the Doppler shift corresponding to the 
electron drift velocity. Furthermore, instabilities will not develop at these short wave- 
lengths even for high electron drift velocities since there is essentially no coupling be- 
tween the electrons and ions over distances much less than the Debye length. 

When the wavelength is much larger than the Debye length {a » 1), the differential 
backs cattering cross section becomes unsymmetric about the frequency of the incident 
wave if the electrons have a net drift velocity. The spectrum develops a sharp peak which 
corresponds to scattering from the progressively more weakly damped ion acoustic waves 
as the drift velocity approaches its critical value. 

The total backscattered power is obtained by numerical integration and is greatly en- 
hanced as the drift velocity approaches the onset of instability. 


Backscatter Cross Section With Electron Drift 

The electrons and ions are both assumed to have Maxwellian velocity distributions, 
with the electrons having a net fluid velocity in the direction of propagation of the 
incident wave. The power spectrum of the components of the stochastic force F , which 
may be thought of as driving the fluctuations, is given by equation (B3) in a reference 
frame in which the s species has zero net fluid velocity. Thus, in the laboratory refer- 
ence frame, the spectrum of F g is Doppler shifted by a frequency kV^ and becomes 

2 

/|F?(ld) F?(ld) |\ do. = k * Ts (zf, - Z? *) did (55) 

X ’ 7 2„L 3 «" - kV d> V 

where Z is the impedance tensor as measured in the laboratory frame. 
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This shows that, for scattering from a plasma in which both species of particles have 
Maxwellian distributions and are drifting with the same velocity, the scattered signal is 
merely Doppler shifted by a frequency equal to kV^. It will be shown that this is also the 
result for the case where the electrons are drifting through the plasma while the ions re- 
main stationary if the wavelength is much smaller than the Debye length. 

The fluid velocity of the electrons is now assumed to be the sum of the constant drift 

1 Cl) f" lk7 

velocity plus a small part proportional to e ; that is 


iwt-ikz 


u_ = Vj + u e 
^6 


(56) 


When the transport equations are linearized, with the zero-order terms assumed to add 
to zero and the impedance tensor Z g solved for in the same manner as done previously, 
the form remains identical with that given by equation (C17) or equation (46) except that 
the frequency c o is replaced by its Doppler shifted value (co - kV^) . Or, equivalently, 

B e is replaced by (6 e - x e ) where x e is the electron drift velocity normalized with re- 
spect to the electron thermal velocity; that is, 



(57) 


The differential backs cattering cross section may formally be written in the same 
form as equation (B14) but with the tensor B now defined as 
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in place of equation (B16). The tensors Z Q and remain of the same form as given by 
equation (C17) except that, in Z g , must be replaced by its Doppler shifted value, 
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The approximate expression for /3 = 90° (eq. (51)) is likewise altered to read 
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( 59 ) 


with z e /P e given by equations (49) and (52) but with fl again to be replaced by 
(®e * e )- 

For the case of a — 0 (Debye length — °°) the differential cross section (eq. (59)) 
reduces to 


a b (co - o>) du> = -Nr 2 ^ — (60) 

b ° W " “ kV d 

which is the same spectrum as for the zero drift case except that the scattered signal is 
Doppler shifted by kV^. This is what one would expect for very large Debye lengths 
since the electrons then move independently of the ions, and the ions have no influence 
over the electrons, which actually cause the scattering. 

As the magnetic field strength B q goes to zero, P g = P. = 1 and the z g are again 
given by equation (27) (with Q & replaced by 6 e - X e in z e )- 

It must be emphasized that in order to legitimately use these expressions for the 
backs cattering cross section, one must be certain that the drift velocity of the electrons 
is less than the critical velocity at which the two-stream instability sets in. 


Numerical Evaluation of Cross Section 

The differential backscattering cross section is evaluated in figure 11 for the case of 
zero net electron drift velocity and for two drift velocities less than the critical drift 
velocity which corresponds to x e = 0. 0055 for this particular set of parameters. This 
figure shows that the spectrum changes from its symmetric shape for X e = 0 to an un- 
symmetric form when the electrons are drifting with respect to the ions. As the drift 
velocity approaches the onset of instability, the spectrum becomes very sharply peaked 
at a frequency corresponding to the ion acoustic wave in the plasma. 

The total scattering cross section is enhanced when the electrons are drifting with 
respect to the ions. As the electron drift velocity approaches the point at which the insta- 
bility commences, the predicted total cross section becomes infinitely large. (In prac- 
tice, the scattered power is limited by nonlinear effects. ) This effect is shown in figure 
12 for four values of the normalized collision frequency i/a . 
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Figure 11. - Normalized Thomson scattering cross section 
as function of Doppler shift for various values of nor- 
malized electron drift velocity, a = 170; normalized ion- 
neutral collision frequency, tfq, 2; normalized electron- 
neutral collision frequency, $ e , 0.2; electron to ion 
temperature ratio, 1 ; normalized ion cyclotron frequency, 
(p- [t 0. 12; angle between incident wave vector and magnetic 
field, p, 90°; C e = D e = 1- Cj = Dj = 2. 

Normalized ion- 



Normalized electron drift velocity, X e 


Figure 12. - Total scattered power as function of normalized 
drift velocity for several values of normalized ion-neutral 
collision frequencies. Magnetic-field strength, zero; 
a= 12.7; normalized electron -neutral collision frequency, 
^ e , 0. 1 i/ty, electron to ion temperature ratio, 10; ion mass, 
mj, 31 atomic mass units; Cj = Dj » 2; C e = D e = 1. 
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SUMMARY OF RESULTS 


The theory of Thomson scattering has been extended to include the effects of colli- 
sions for scattering from a plasma with unequal electron and ion temperatures including 
the influence of a constant magnetic field and electron drift. The theory presented here 
is based on the two -fluid continuum equations for a plasma and should provide a good de- 
scription when the wavelength of the probing wave is larger than the ion-neutral mean 
free path. 

When the ions and electrons are not at the same temperature, the total backscattered 
power is shown to be dependent on the ion-neutral and electron-neutral collision frequen- 
cies, as well as the magnitude and relative angle of the magnetic field if one is present. 
For no magnetic field and electron to ion temperature ratio T e /T. greater than or equal 
to 10, the total backscattered power is shown to be substantially less than that of colli- 
sionless calculations but always greater than a cross section that varies with T_/T. as 
- 1 e i 

(1 + T /T.) . For example, with o » 1 {a = X./47TX,-, where X^ is the wavelength of 

“ 1 o u o 

the incident wave and is the Debye length) T g /T^ = 20, and the ion-neutral mean free 
path X. ~ 1/5 X , the normalized total backscattered power is about 1/7 of that predicted 

111 ° i i 

for no collisions and about times the value given by (1 + T /T-) . 

2 e 1 

As the electron -neutral and ion-neutral collision frequencies v and v. increase, 

en in 

the scattered power approaches the (1 + T /Tj)“ curve, although it does not reach it 
even in the limit v , v ■ — Also, for T /T- > 5, a change in the ratio of the colli- 

sion frequencies v / f- n can have a large effect on the total scattered power although 
the shape of the spectrum is essentially unchanged. 

A constant magnetic field affects Thomson scattering only when the field is normal 
to, or nearly normal to, the direction of propagation of the incident wave. For T/T- 
in the common experimental range of 1 to 10, the total scattered power is increased by 
the magnetic field. In fact, the scattered power may increase with T/T. rather than 
falling off in accord with the usual (1 + T /T-)“ relation. For instance, with both the 
ion gyroradius and the ion-neutral mean free path equal to 1/5 X , with the angle between 
the incident wave and the magnetic field ft equal to 90°, and with T e /T^ ~ 5, the total 
backscattered power is almost five times that calculated for the field-free case. Increas- 
ing the collision frequencies lessens the influence of the field and causes the scattered 
power to become closer to the field -free case. In general, the total scattered power is a 
function of the collision frequencies and the magnetic field strength except for the special 
case of T e = T- when it is totally independent of both the collisions and the magnetic 
field. When electrons are drifting with respect to the ions, only the component of the 
drift velocity in the direction of the incident wave affects the backscattered signal. The 
spectrum is no longer symmetric about the incident frequency (as it is for no drift). If 
the Debye length is large compared with the wavelength of the incident wave, the only 
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effect on the spectrum is a Doppler shift corresponding to the component of the drift 
velocity in the direction of the incident wave. However, if the Debye length is small com- 
pared with the wavelength, an instability of the ion acoustic wave develops as the drift 
velocity reaches a certain critical value. 

As the drift velocity increases, the spectrum of the backscattered signal develops a 
sharp peak at a frequency shift (from the incident wave frequency) corresponding to the 
weakly damped ion acoustic wave. The total scattered power is enhanced by a nonzero 
drift velocity and, as the electron drift nears the critical velocity, becomes very much 
greater than that calculated for a quiescent plasma. For example, with no magnetic 
field, Tg/T. = 10, a » 1, and i X q , an electron drift velocity equal to 1/20 of the 

electron thermal velocity causes the total backscattered power to be increased by a factor 
of 14 over the zero -drift case. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, February 26, 1971, 

120-26. 
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APPENDIX A 


SYMBOLS 



C 


c s 

D 

D s 

E 

e 


ft;S 

H 

n 

i 

Im 

i 


J 



k 


k 


o 


k 


s 


L 


m 


s 


N o 


see eq. (50) 
see eq. (C13) 

constant magnetic flux density 

see eq. (B16) 

constant of order 1 or 2 

velocity of light 

see eq. (C9) 

electric flux density 

constant of order 1 or 2 

electric field intensity 

electronic charge 

applied force 

see eq. (C12) 

magnetic -field intensity 

Planck's constant/2-?r 

generalized system response 

imaginary part 

>/-l 

current density 

see eqs. (40a) and (40b) 

wave vector of scattering wave 

wave vector of incident wave 

wave vector of scattered wave 

linear dimension of scattering volume 

mass 

ambient electron number density 
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N s 

n 

P„ 


Q 

% 

Re 


os 


U 


~s 


V d 

W 


h’l'h 

z s 

a 

P 

P' 

Es 

r 

A s 

AN 


number density 

index of refraction (for scattering wave) 
see eq. (52) 

traceless pressure tensor 

scalar pressure 

see eq. (39) 

charge 

real part 

classical electron radius 
see eq. (CIO) 
equilibrium temperature 
temperature 
see eq. (Cll) 
plasma sound speed 
mean velocity 
generalized force 
electron drift velocity 
energy in L 
see eq. (47) 
iwe 0 /e 2 

generalized admittance tensor 
normalized impedance tensors (see eq. (Bll)) 
normalized impedance function (see eq. (27)) 
l/(kX D ) 

angle between k and 13 
see eq. (35b) 
flux density 

square root of electron to ion mass ratio 
see eq. (C19) 

fluctuating part of electron number density 
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Ae 


6 ij 


e(r,t) 

e„ 


V, 

e 

e i 

K 

X 


OS 


D 


m 


os 


M 

^sn 


°b<"o - "> 
°tot 


% 

oo 

oo 

e 


fluctuating part of permittivity 
Kronecker delta 
permittivity of medium 
permittivity of free space 
coefficient of viscosity 
scattering angle 

normalized Doppler shift frequency (see eq. (29a)) 

Boltzmann's constant 

electron Debye length 

ion -neutral mean free path 

thermal conductivity (with B c = 0) 

wavelength of incident wave 

thermal conductivity (B ^ 0) 

electron to ion temperature ratio 

effective collision frequency for momentum transfer 

effective collision frequency for energy transfer 

angle between incident electric field vector and k 

rx/ S 

mass density 
see eq. (C20) 

differential backs cattering cross section 
total backs cattering cross section 
normalized cyclotron frequency (see eq. (29b)) 
normalized electron drift velocity (see eq. (57)) 
normalized collision frequency (see eq. (C7) 
cyclotron frequency 
Doppler shift frequency 
plasma frequency 
frequency of incident wave 
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Subscripts: 

e 

i 

s 

so 

n 

i,j,k,l,m,n 

Superscripts: 

* 

t 


electrons 

ions 

either electrons or ions 

zero -order quantity; either electrons or ions 
neutrals 

components of vectors or tensors 

complex conjugate 
Hermitian conjugate 
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APPENDIX B 


BACKSCATTERING CROSS SECTION FOR ANISOTROPIC MEDIUM 


A general tensor admittance function Y„ is defined in the same manner as 
Dougherty and Farley (ref. 3): 


r a = Y 0 

~ S RiS 


F s 


(Bl) 


where F is an applied force and r_ is the resulting flux density. The rate of change 

~ S o ^ S 

of energy in a volume L° is 


— = L 3 r . F (B2) 

dt ~ s ~ s 


Applying the fluctuation -dissipation theorem (ref. 17) to the s species (for kT » hu>) 
gives the cross spectral density of the components of F 

<-o S 

<|Ff MFjV) |) (to - flsL to (B3) 

2itIj ^ 


This may be thought of as the spectrum of the stochastic forces that are driving the s 
species when no interaction due to the electric fields is taken into account. The cross 
spectral density is the Fourier transform of the cross correlation function and is defined 
by 



e iwT do> = (Ff*(t)F*(t + rj> 


(B4) 


With the inclusion of the force due to the electric field, the equation of motion for 
the electrons (and ions) is 


= Js 


< F s + 




(B5) 
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Maxwell’s equations provide another relation 


with 


r 


e 



• (eE) 



1 - n 2 0 0 

0 1 - n 2 0 

. 0 0 1 . 


(B6) 


(B7) 


and where the coordinate system has been chosen so that the wave vector k lies along 
the z axis . 

Using equations (B5) and (B6) and assuming that the stochastic forces F e and f\ 
are uncorrelated lead to the cross spectral density of the electron fluctuations 


where 



d - = Y ij Y 


* 

fm 




Y 


» * 
mn 



(B8) 


Y = Y q • (Y a + Y. + Y t ) _1 • (Y. + Y t ) 
Y’ = (y. + Y t ) _1 • Y. 


(B9) 

(BIO) 


Normalized impedance tensors Z , Z and Z T are defined by 

Rd S ftii 

iojN 


Y = — 

2 


z; 1 

PiiS 


k*xT c 


iwN 
Y = -Z' 


k 2 K T 


(Bll) 


icuN 


Y i = 


k 2 /cT 


° Z T _1 
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The differential backs cattering cross section is proportional to the power spectrum 
of the Fourier component of the electron number density fluctuations with wave vector 
2k (ref. 3). 

rsjVJ 


CT b (u> o - w) dco = L 3 r 2 ^AN(2k Q , w) | 2 ^ dw (B12) 

The backs cattering cross section is found by using equations (B7) to (B12) together 
with the continuity equation for the electrons 


3AN 

at 


+ v • 



(B13) 


to be 


where 


and 


a b (u>0 - co) do) = 


Nr 2 
o e 

27TiOJ 




(B14) 


Z 


= Z g + 


(mi r 1 - 


Zf 1 ) 


lx-1 


(B15) 


B 




zf 1 )" 1 


(B16) 
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APPENDIX C 


CALCULATION OF IMPEDANCE TENSORS 

In this appendix the normalized impedance tensors Z g and Zj are calculated for 
the electrons and ions considered separately; that is, without including the force due to 
the self-consistent electric field. The impedance tensor Z a for the s species is de- 
fined by 


icuN 

z o • r = — f 


k 2 *T 


s 


(Cl) 


where r is the flux density resulting from an applied force F_ that is proportional to 

iu>t-ikz~ s ~ s 

e 

The approach of Tanenbaum (ref. 11) is followed where the transport equations for 
mass, momentum, and energy are used. For small amplitude disturbances in a weakly 
ionized gas with the mean velocity of the neutral particles u equal to 0, it is found to 

r<jTl 

first order in the fluctuating quantities, that 


d _Ps 

at 


+ • u^, = 0 

^so ~ s 


(C2) 


3u Vp F q V • P^ 

(a s XB 0 ) + — 
dt P so m s m s Ac 


sn~s 


(C3) 


'so 



at 


5^ p so s 

3 9t 
^so 


-LV 2 T = -N KV ' (T - T J 
2 s s o sn ' s so 7 


(C4) 


where p , q g , m g , p g , and T g are the mass density, charge, mass, pressure, and 

temperature, respectively, for species s; B Q is a constant magnetic field; and v gn 

and v' = 2m v „/(m + m ) are, respectively, the effective collision frequencies for 

momentum and energy transfer with the neutrals. The form of the traceless pressure 

tensor and heat flow term developed by Chapman and Cowling (ref. 24) is used. For all 

i c±)t ik7 

fluctuating quantities varying as e , the divergence of the traceless pressure ten- 

sor P g may be written in the form 



Normalized frequencies for the Doppler shift frequency, collision frequency, and 
cyclotron frequency are defined by 


e s = 


to / m s 


1/2 


k \2 K T 


(C6) 


K = 


v sn / m s 


1/2 


k \2/cT. 


(C7) 


% m s 


1/2 


k \2*T 


(C8) 


Then T g is a matrix whose elements are functions only of 






(C9) 


s 


s 


= + 


'* s \ . a 

_ ) sin /3 

vV 


(CIO) 


where /3 is the angle between the incident wave vector k and the magnetic field B. . 

'"“•'O 

The upper sign is to be taken for the electrons, and the lower sign for the ions. The ten- 
sor is then given by 



(Cll) 


where 
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G„ = 


1 + 



1 +4c' 


2s, 


1 + 4c‘ 



4s' 


1 + 4c' 


2s' 


1 + 


1 +4c' 


2s, 



(C12) 


and 


A„ = 


1 0 


0 1 


0 0 


4 

3 

2 

3 


c s 
s s 


1 +4ct 


1 +4ct 


4 

3 


(C13) 


The thermal conductivity is given by 


1 + ( tl \ cos 2 /3 


X s = X os 


v^= 




(C14) 


The coefficients for viscosity and thermal conductivity for charged particles in a weakly 
ionized gas which appear in equations (C5) and (C14) are given by 


08 D v 

s sn 


(C15) 


and 


\ = 
os 


15/cp c 


(C 16) 


4m t ,C c , v on 
s s sn 
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where D g and C g are constants (of order 1 or 2) that depend on the interparticle force 
law and the mass ratio between the electrons or ions and neutrals . 

By using the mass and energy transport equations, we may write the momentum 
equations in the same form as equation (Cl) with the impedance tensors Z given by 


ie s 

z a = — 5 


D xly 
s*s 


T + 


+210 0 cos p 


+2i0 s 0 s cos ^ 


s s 
h s 


+2i0 g e s sin /3 


±2i0 s 0 s sin ^ 


A s + h s 


(C 17) 


where 


and 


h s = 2i ^s - 2e s 


(C18) 


l + if^ 

\3cr„ . 


A s = 


1 + if" 


(C 19) 


5k 2 /cT 


CT s = v sn + - 


1 +(— \ cos^jS 


2< Vsn m s 


1 




s/ J 


(C20) 


Now equation (C17) may be used to evaluate Z g and Z^ in terms of plasma parameters. 
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